N— ( 


CO 

to 

i> 

O: 


EXPERIMENTAL  PULSED  LASER  REMOTE  CROSSWIND 
MEASUREMENT  SYSTEM  - 
FEASIBILITY  STUDY  AND  DESIGN  (PART  V) 


J.  FRED  HOLMES 
J.  RICHARD  KERR 
RICHARD  A.  ELLIOTT 
MYUNG  H.  LEE 
PHILIP  A.  PINCUS 
MICHAEL  E.  FOSSEY 


DEPARTMENT  OF  APPLIED  PHYSICS  &  ELECTRONIC  SCIENCE  D  D  C 
THE  OREGON  GRADUATE  CENTER  fTMTirTnrpnn  P 

19600  N.W.  WALKER  ROAD  !  H  |l  Li  P-M  rJlUilU- 

BEAVERTON,  OREGON  97005 


DEC  5  3ST9 


SEPTEMBER  1978 


iitisiEirDnE 

^  A 


s 


AKiiV 


ARMAMLNI  RESEARCH  AND  DEVELOPMENT  COMMAND 
FIRE  CONTROL  AND  SMALL  CALIBER 
WEAPON  SYSTEMS  LABORATORY 
DOVER,  NEW  JERSEY 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED., 


10 


,  } 

'J 


r 


v 


NOTICES 


Disclaimers 
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INTRODUCTION 


In  FY  74  Frankford  Arsenal  through  the  Atmospheric  Sciences 
Laboratory  (ECOM)  White  Sands  Missile  Range,  New  Mexico,  initiated  a 
program  to  develop  a  remote  crosswind  measurement  system  for  ballistic 
applications.  A  remote  crosswind  system  is  needed  that  will  measure 
the  average  crosswind  between  gun  and  target  and  automatically  feed 
this  data  to  an  onboard  ballistics  computer  for  fire  control  correc¬ 
tion  of  the  crosswind  error.  This  error  is  quite  significant,  repre¬ 
senting  one  of  the  largesl  remaining  sources  of  combat  vehicle  fire 
control  error,  especially  with  medium  and  slow  muzzle  velocity 
ammunitions  (ref.  1) . 

In  an  effort  to  meet  this  need,  a  pulsed  laser,  remote  crosswind 
measurement  system  is  being  developed  at  the  Oregon  Graduate  Center. 
The  present  system,  although  not  optimized  and  fully  tested,  appears 
to  provide  an  accurate  indication  of  the  path  averaged  crosswind  be¬ 
tween  the  transmitter  and  target.  It  still  must  be  demonstrated  that 
the  system  will  work  accurately  at  longer  ranges  and  at  low  turbulence 
levels  and  the  limits  of  accurate  system  operation  must  be  defined  and 
extended.  Work  on  these  tasks  and  others  will  be  accomplished  during 
the  next  contract  period.  The  tasks  include: 

Measurement  Pulsed  System 


Non-Time  Delayed  Statistics  in  Extended  Ranges 
Time  Delayed  Statistics 

Wind  Measurement,  Variety  of  Atmospheric  Conditions 

Equipment  and  Software  Improvement 

Wind  Calculation  Between  Pulse  Pairs 
Improved  Diode  Arrays 
Reduction  of  Pulse  Separation 

Self-adaptive  to  Signal  Strength,  Range  and  Wind 
Error  Reduction  Through  Preprocessing  Data  and  Improved 
Wind  Algorithms 

Beam  Alignment  Effects 

Analytic  Formulation 
Definition  of  Laser  Requirements 

Coherence  Effects 


Analytic  Formulation  of  Partially  Coherent  Source  Problem 
Definition  of  Laser  Requirements 


1 


System  Optimization 


Error  Analysis 
Numerical  Studies 


BACKGROUND  AND  SUMMARY 

The  feasibility  of  remotely  determining  the  transverse  wind  vel¬ 
ocity  with  a  point-to-point,  cw  laser  system  is  well  established  (refs. 
2-4).  Tn  more  recent  work  (refs.  5-1.1),  investigators  considered  the 
extension  the  technique  to  the  use  of  multiple  cw  lasers  and  folded 
paths  as  ue±'  as  ambient  light  and  the  use  of  techniques  to  overcome 
saturation.  1,  '  technique  basically  consists  of  relating  the  scin¬ 
tillation  patten,  on  the  receiver  to  the  frozen-in,  transverse  motion 
of  the  turbulence  structure  along  the  propagation  path.  The  most  suc¬ 
cessful  cw  method  that  ''as  been  used  involves  the  measurement  of  the 
slope  of  the  log  amplitue  covariance  function  at  zero  time  delay  (ref. 
2)  . 


Using  active  cw  sources,  an  ’cceptable  signal  to  noise  ratio  can 
be  achieved  at  reasonable  transmitter  power  levels.  However,  the  pro¬ 
posed  application  of  the  work  sponsor*. 1  under  this  project  requires  a 
single-ended  system  operating  against  a  on-cooperative  target  which 
can  be  viewed  as  a  perfectly  diffuse  surfa^'»  with  reflectivity  of 
approximately  0.1  (ref.  12).  Under  this  cons  Taint,  an  active  cw 
system  would  require  excessive  transmitter  powe,  levels. 

The  primary  noise  sources  in  the  system  are  the  poisson  fluctu¬ 
ations  of  a  detector  current  and  thermal  noise  in  the  v’ectronics. 

For  a  given  pulse  energy,  as  the  pulse  width  decreases  tu *  instantan¬ 
eous  power  increases.  Consequently,  che  detector  signal  cu.  'rent  in¬ 
creases  and  eventually  its  Poisson  fluctuation  becomes  the  dou 'nant 
noise  source.  If  the  receiver  is  gated  synchronously  with  the  i.  -*ns- 
tnitteu  pulses,  it  can  be  made  signal  shot  noise  limited.  Therefor.,  in 
a  direct  detection  optical  system,  the  signal  to  noise  ratio  for  a 
given  total  energy  increases  as  the  pulse  duration  is  made  shorter 
(ref.  13)  until  the  signal  shot  noise  limit  is  reached.  Consequently, 
a  pulsed  system  can  be  made  considerably  more  sensitive  than  a  cw 
system. 

In  order  to  employ  a  pulsed  laser  and  still  make  a  temporal 
measurement  as  required  for  wind  velocity,  a  double-pulse  technique 
has  been  used  (refs.  14,  15,  16)  whereby  two  Q-switched  laser  pulses 
are  generated  at  a  time  interval  on  the  order  of  a  millisecond.  This 
interval  appears  to  be  appropriate  for  the  employment  of  the  covariance 
"slope  at  zero  time  lag"  technique  (ref.  2)  and  is  compatible  with 
current  double-pulse  laser  technology.  Initially,  a  ruby  laser  was 
used  as  a  pulsed  source. 
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Since  the  pulsed  technique  utilizing  a  ruby  laser  precluded  time 
averaging  of  the  statistical  quantities  measured,  spatial  averaging 
had  to  be  employed.  Therefore,  an  array  of  detectors  was  used  in  order 
to  generate  a  sufficient  number  of  samples  of  the  instantaneous  scin¬ 
tillation  pattern  to  use  spatial  averaging.  To  the  extent  that  the 
statistics  of  this  pattern  are  invariant  with  lateral  translation  at 
the  receiver,  such  spatial  averaging  will  give  the  desired  result;  and 
the  pulsed  technique  is  fundamentally  similar  to  the  cw,  time-averaging 
approach. 

Initially  the  work  on  the  pulsed  laser,  remote  crossvind  measure¬ 
ment  system  was  based  on  a  technique  and  theory  developed  for  a  cw 
system  utilizing  a  laser  directly  as  a  source.  This  technique  and 
theory,  which  is  based  on  log-amplitude  statistics,  has  been  verified 
experimentally  and  is  valid  for  the  cw  configuration.  The  basic  cw 
technique  is  applicable  to  the  pulsed  system;  however,  the  detailed 
theory  does  not  apply  to  the  pulsed  configuration.  This  is  not  the 
result  of  a  difference  in  the  propagation  characteristics  between  the 
pulsed  and  cw  sources,  but  is  caused  by  the  physical  differences  in 
the  two  systems. 

The  source  for  the  pulsed  system  is  a  laser  illuminated  incoher¬ 
ent  target  yielding  a  return  path  source  with  an  undefined  size  which 
at  least  close  to  the  target  gives  rise  to  a  temporally  coherent, 
spatially  incoherent,  Rayleigh  distributed  field  amplitude.  This  is 
in  contrast  to  the  point  source  used  in  the  earlier  active  cw  systems 
and  the  temporally  incoherent  cw  sources  used  in  later  active  and  pas¬ 
sive  systems  which  give  rise  to  log-normal  amplitude  statistics. 

Consequently,  the  experimental  work  with  the  ruby  transmitter  con¬ 
centrated  largely  on  tests  where  the  size  of  the  source  (at  the  target) 
could  be  controlled  (ref.  16).  This  allowed  a  better  comparison  with 
available  theory  by  decoupling  from  the  experiment  the  effects  of  the 
atmosphere  on  the  laser  beam  as  it  propagates  from  the  transmitter  to 
the  target.  The  controlled  source  size  was  accomplished  by  utilizing 
small  disks  of  scotchlite  on  a  black  background  as  a  target. 

The  results  of  these  tests  indicated  the  need  both  for  a  better 
pulsed  source  and  for  a  theory  applicable  to  the  pulsed  system. 

A  ruby  laser  was  originally  chosen  as  a  source  because  of  cost 
considerations  and  the  realization  that  energy  levels  of  the  order  of 
a  joule  or  greater  might  be  required.  Subsequent  analysis  and  ex¬ 
perimental  results  indicated  that  energy  levels  of  the  order  of  100 
millijoules  would  yield  an  adequate  signal  to  noise  ratio  in  the  sys¬ 
tem.  Consequently,  the  system  was  redesigned  to  utilize  a  Nd:YAG 
laser  operating  at  1.06  microns.  It  has  the  following  advantages  over 
the  ruby  system:  much  smaller;  more  rugged;  less  sensitive  to  tempera¬ 
ture  changes;  lighter;  less  power  required;  much  less  Radio  Frequency 
Interference  (RFI) ;  convection  cooled;  space-time  averaging  possible. 
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Most  of  the  advantages  of  the  NdrYAG  system  accrue  from  the  much 
lower  threshold  of  the  lasing  medium.  Consequently,  it  requires  a  much 
lower  input  energy  per  pulse  (on  the  order  of  10  in  lieu  of  1,000 
joules) .  This  allows  the  YnG  laser  to  be  double  pulsed  at  a  much  high¬ 
er  rate  than  the  ruby  laser  and  makes  space-time  averaging  possible. 

In  addition,  the  lower  input  energy  greatly  reduces  RFI  between  the 
laser  and  the  receiver. 

The  use  of  space-time  averaging  allowed  the  size  of  the  optical 
receiver  array  (spatial  averaging)  to  be  reduced  from  64  elements 
(8  *  8)  to  16  elements  (2*8).  The  laser  can  produce  a  double  pulse 
once  every  100  milliseconds ,  and  the  number  of  double  pulses  used  for 
a  crosswind  measurement  can  be  adjusted  from  one  to  forty.  An  addi¬ 
tional  advantage  of  the  time  averaging  is  that  it  tends  to  eliminate 
any  intermittency  of  turbulence  problems.  The  YAG  system  was  construc¬ 
ted  and  at  the  start  of  the  present  report  period  was  ready  for  verifi¬ 
cation  tests  and  debugging.  This  phase  of  the  development  required  al¬ 
most  10  months  and  consisted  of  system  problem  identification  and 
solution.  Some  of  the  problems  encountered  and  solved  were:  the  need 
to  directly  measure  the  average  received  intensity  as  opposed  to 
calculating  it;  the  need  for  an  effective  method  of  aligning  the  two 
laser  beams  in  the  transmitter  so  that  they  are  colinear;  the  need  for 
a  method  of  normalization  to  compensate  for  turbulence  strength  and 
laser  coherence;  noise  and  crosstalk  in  the  receiver;  the  need  for  an 
adjustable  time  delay  between  pulses;  and  the  need  for  accurate,  auto¬ 
matic  gain  calibration  of  the  16  receiver  channels.  The  pulsed  wind 
sensor  is  now  fully  operational  and  is  providing  significant  data 
with  which  to  evaluate  and  optimize  the  overall  system. 

Since  a  theory  applicable  to  the  pulsed  system  was  not  available, 
a  new  theoretical  formulation  for  the  statistics  of  the  received  in¬ 
tensity  based  on  the  extended  Huygens -Fresnel  principle  was  developed 
(refs.  16-18).  Our  first  theoretical  effort  invoked  the  assumption 
of  Gaussian  field  statistics  at  the  receiver.  Because  of  this  assump¬ 
tion,  only  phase  perturbation  effects,  which  show  up  as  a  dependence 
on  the  transverse  correlation  scale  (p  ),  v?hich  appear  in  the  result. 
This  was  initially  thought  to  be  a  goo3  assumption,  at  least  for  low 
turbulence  levels,  on  the  basis  that  the  relation  between  the  two  tur¬ 
bulence  scale  sizes  (p  and  /L/k)  was  expected  to  be  multiplicative 
and,  therefore,  the  smaller  scale  -  which  is  usually  p  in  the  present 
application  -  would  be  dominant. 

Our  cw  experiments  (ref.  17)  confirmed  that  at  low  turbulence 
levels,  the  Gaussian  assumption  is  indeed  correct  for  the  variance  and 
covariance.  However ,  at  intermediate  turbulence  levels  the  scale 
size  approaches  /L/k.  This  suggested  that  the  log-amplitude  term, 
which  is  responsible  for  this  scale  size,  plays  an  important  part  in 
the  statistics  of  the  received  intensity  and  particularly  with  regard 
to  remote  wind  sensing. 
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Consequently,  a  complete  theoretical  formulation  including  the  log- 
amplitude  term  was  developed  (ref.  17).  It  shows  that  the  relationship 
between  the  two  scale  sizes  is  additive  rather  than  multiplicative  as 
indicated  by  the  cw  experimental  data.  This  means  that  >'L/k  scale  size 
can  be  dominant  for  our  application  and  that  the  same  mechanism  which 
has  been  used  as  a  basis  for  the  successful  cw  wind  sensing  work  per¬ 
formed  by  NOAA  is  available  for  utilization  by  the  pulsed,  diffuse- 
target  case. 

During  the  present  contract  period,  additional  data  was  obtained 
from  our  cw  experiments  to  increase  our  confidence  in  the  complete 
theoretical  formulations.  In  addition,  the  cw  experimental  site  was 
relocated  adjacent  to  the  pulsed  experimental  facility.  This  required 
a  major  effort  including  constructing  a  new  building,  and  moving  and 
realigning  the  cw  equipment.  The  reason  for  this  action  was  to  pro¬ 
vide  a  propagation  path  for  the  cw  experiment  where  the  wind  is  pri¬ 
marily  across  the  path.  This  is  needed  in  order  to  measure  the  time- 
delayed  statistics  for  comparison  with  the  theory.  These  measurements 
have  been  completed  and  good  qualitative  agreement  with  the  theory  was 
obtained. 

Analytic  work  during  the  present  contract  period  was  concentrated 
on  numerical  work  in  support  of  the  pulsed  and  the  cw  experimental 
efforts.  Specifically,  more  efficient  techniques  were  developed  for 
obtaining  numbers  from  the  theoretical  formulations  which  included 
algorithms  for  evaluating  the  two  source  wave  structure  function  and 
two  source  log-amplitude  covariance  function. 

SYSTEM  CONSIDERATION 

Several  new  techniques  were  developed  for  improving  system  per¬ 
formance  during  the  contract  period.  These  included:  finding  a  method 
of  aligning  the  laser  beams;  adding  hardware  and  software  for  direct 
average  measurements;  introducing  a  time  delay  between  when  the  laser 
starts  pulsing  and  when  data  is  taken;  and  gain  leveling. 

The  two  laser  beams  used  in  the  transmitter  must  be  colinear  in 
order  to  measure  the  time  delayed  statistics  and  hence  the  crosswind. 
Initially,  the  beams  were  aligned  by  observing  them  with  an  infrared 
viewer  at  the  target  and  adjusting  the  transmitter  until  the  spots 
overlap.  This  was  difficult  and  highly  inconvenient  because  the  target 
spots  could  be  resolved  only  during  very  low  turbulence  and  ambient 
light  conditions,  and  the  observer  and  adjuster  were  physically  sep¬ 
arated  by  500  meters. 

In  order  to  overcome  this  problem,  a  method  of  checking  and  ad¬ 
justing  beam  alignment  within  the  confines  of  the  field  site  building 
was  developed.  The  method  involves  frequency  doubling  the  1.06  micron 
radiation  from  the  transmitter  to  obtain  visible  light  (green)  and 
then  projecting  the  far  field  pattern  on  a  screen.  The  arrangement  is 
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illustrated  in  figure  1, 

A  secondary  benefit  of  converting  to  visible  light  for  alignment 
was  the  discovery  that  during  the  first  few  seconds  of  operation,  the 
laser  beam  is  highly  distorted.  We  are  now  compensating  for  this  by 
using  a  longer  burst  of  pulses  and  ignoring  the  data  for  the  first  2 
seconds.  A  noticeable  improvement  was  achieved. 

Because  the  laser  used  in  the  system  is  temporally  quite  incoher¬ 
ent,  the  atmospherically  induced  fluctuations  in  received  intensity 
which  carry  the  wind  information  are  only  a  fraction  of  the  mean  value 
of  received  intensity.  Consequently,  the  mean  value  must  be  measured 
with  a  high  degree  of  precision  in  order  not  to  create  a  large  error 
in  the  crosswind  measurement.  This  requires  that  an  aperture  averaged 
measurement  of  the  mean  intensity  be  made  and  that  careful  attention 
be  paid  to  gain-leveling  the  sixteen  receiver  channels. 

We  had  anticipated  that  averaging  the  sixteen  receiver  channels 
might  not  yield  a  mean  intensity  measurement  of  sufficient  accuracy; 
and  consequently  had  made  provision  for  adding  a  beam  splitter  and 
mean  value  receiver  in  the  optical  train.  Our  subsequent  experimental 
work  confirmed  the  need  for  a  direct  average  measurement.  However, 
when  we  tried  to  implement  the  beam  splitter  technique,  our  efforts 
were  frustrated  by  the  fact  that  the  laser  had  been  designed  such  that 
the  two  pulses  were  orthogonally  polarized. 

With  a  perfectly  diffuse  target,  the  return  should  be  randomly 
polarized  and  independent  of  the  polarization  of  the  incident  beam. 

In  practice,  it  was  found  that  real,  diffuse  like  targets  retain  some 
of  the  polarization  of  the  incident  beam.  This  coupled  with  a  polari¬ 
zation  sensitive  beam  splitter,  convinced  us  to  terminate  the  beam 
splitter  approach  and  to  build  a  small  separate  receiver  to  monitor 
the  mean  intensity  from  each  pulse.  In  a  field  version  of  the  system, 
the  beam  splitter  approach  should  be  used  to  save  space.  Consequently, 
the  next  version  of  the  laser  transmitter  should  either  be  circularly 
polarized  or  have  both  pulses  linearly  polarized  in  the  same  direction. 

Gain  leveling  of  the  receivers  is  done  automatically  by  the 
processor.  In  order  to  gain  level  the  receivers,  a  uniform  intensity 
is  required  to  be  incident  on  the  receiver  optics .  This  is  accom¬ 
plished  by  illuminating  a  target  some  distance  from  the  receiver  with 
the  laser  transmitter.  The  received  intensity  for  any  particular  pulse 
will  vary  from  detector  to  detector  due  to  turbulence.  However,  on  the 
average  it  will  be  the  same  value.  Consequently,  by  taking  many 
samples  (600  presently)  and  averaging,  the  relative  gain  of  each  of 
the  sixteen  channels  and  the  average  receiver  is  determined.  Separate 
gain  constants  for  each  of  the  two  pulses  in  a  pulse  pair  are  deter¬ 
mined  and  stored  by  the  computer.  This  procedure  yields  a  gain 
leveling  good  to  about  one  percent  which  appears  to  be  adequate  for 
the  application.  Best  results  are  obtained  when  one  of  the  center 
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Green  Light  (two  beams) 


Figure  1.  Beam  alignment  setup. 


most  elements  in  the  receiver  array  are  used  as  a  reference. 

LASER  BEAM  ALIGNMENT 

In  order  to  measure  the  time  delayed  covariance  from  which  the 
crosswind  is  derived,  two  laser  pulses  separated  in  time  must  be  used. 
However,  if  the  laser  beams  are  misaligned,  they  do  not  travel  through 
the  same  portion  of  the  atmosphere.  Consequently,  the  atmospherically 
induced  fluctuations  on  the  received  intensity  from  the  two  pulses 
separated  in  time  that  make  up  a  measurement  pair  become  decorrelated . 
Since  it  is  the  correlated  fluctuations  that  contain  the  wind  infor¬ 
mation,  misalignment  of  the  beams  degrades  and/or  prevents  proper 
system  operation. 

A  good  measure  of  the  degree  of  correlation  can  be  obtained  by 
observing  the  autocovariance  function  defined  by 

<(l(r2,t2)  -  l) (l i » t i )  "  l)> 


r  =  0 

where  t2  -  tj  =  T  =  time  delay  between  pulses 

"r  =  ~r2  ~  ri  =  separation  between  detectors 
and  Oj2  is  the  variance 

This  quantity  can  be  thought  of  as  being  proportional  to  the  signal  we 
are  trying  to  measure.  Consequently,  if  it  is  reduced  by  beam  mis¬ 
alignment,  this  is  equivalent  to  reducing  the  signal  to  noise  ratio  of 
the  system  and  thereby  increasing  the  wind  measurement  error.  In 
addition,  as  the  beams  separate,  the  system  calibration  constant 
changes.  This  introduces  additional  measurement  error. 

Figure  2  shows  the  results  of  some  measurements  at  1582  meters. 
The  data  were  obtained  by  misaligning  the  beams  by  a  set  amount  and 
then  measuring  the  autocovariance.  As  predicted,  it  does  fall  off 
significantly  as  the  beams  separate.  If  a  drop  of  10%  in  the  auto¬ 
covariance  can  be  tolerated,  then  the  data  indicate  that  the  beam 
misalignment  on  a  pulse  to  pulse  basis  must  be  no  worse  than  the  laser 
beam  divergence/ 3. 

The  actual  requirement  probably  depends  on  range,  transmitter 
size  and  turbulence  level.  The  beams  must  physically  overlap  to  some 
as  yet  unknown  degree  in  order  for  the  fluctuations  to  be  correlated 
sufficiently  by  the  atmosphere.  Qualitatively  then  it  would  seem  that 
with  the  atmospheric  fluctuations  dominant  (as  opposed  to  speckle), 
there  are  two  regimes  with  respect  to  the  effect  of  range  and  turbu¬ 
lence  on  the  error  caused  by  beam  misalignment: 


C  (r , t) 
N 


8 


Normalized  Auto  Covariance 


Figure  2.  Autocovariance  versus  beam  separation. 
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The  beam  spread  due  to  laser  coherence  and  the  turbulence  is 
(1)  greater  than  the  beam  misalignment  or  (2)  less  than  the 
beam  misalignment.  In  the  first  case  the  sensitivity  of  the 
system  and  the  accuracy  should  increase  with  turbulence  level 
and  range  and  in  the  latter  case  decrease. 

Our  operational  experience  to  date  has  fallen  into  the  second  regime; 
however,  it  would  be  much  more  desirable  to  operate  in  the  first 
regime.  Changes  to  the  laser  transmitter  would  make  this  possible. 

In  the  second  regime,  due  to  the  dependence  of  the  calibration 
constant  on  beam  alignment,  the  fluctuations  in  indicated  crosswind 
should  increase  with  range  thereby  giving  an  increasing  RMS  error  with 
range.  This  has  been  our  experience  to  date  and  is  the  opposite  of 
what  was  expected  prior  to  encountering  the  beam  alignment  problem. 

With  perfectly  aligned  beams,  it  is  expected  that  the  error  will  de¬ 
crease  with  range  (neglecting  S/N  limitations)  due  to  increased  path 
averaging  and  increased  integrated  path  turbulence.  There  may  be  a 
limit  as  to  how  far  in  range  this  argument  might  be  valid  due  to 
saturation  effects;  however,  we  have  not  yet  encountered  anything  ex- 
perientially  that  could  be  reasonably  attributed  to  saturation. 

Another  deleterious  effect  of  beam  misalignment  is  that  it  in¬ 
creases  the  lower  limit  of  turbulence  level  at  which  the  system  will 
operate.  As  the  turbulence  level  decreases,  a  point  will  be  reached 
at  which  the  speckle  field  will  start  to  dominate  (see  the  section  on 
laser  coherence)  and  the  system  will  lose  its  sensitivity  to  the  cross- 
wind.  If  the  beams  are  misaligned,  this  will  occur  prematurely  since 
the  misalignment  desensitizes  the  system  to  the  atmospherically  in¬ 
duced  fluctuations  that  carry  the  wind  information. 

Also,  it  appears  that  a  horizontal  misalignment  of  the  beams  will 
cause  a  constant  bias  to  the  indicated  wind.  The  system  responds  to 
the  shifting  of  the  scintillation  pattern  with  the  wind.  If  the  beam 
for  the  second  pulse  is  displaced  horizontally  from  the  beam  for  the 
first  pulse,  then  the  pat  Lem  will  appear  to  the  system  to  have  shifted 
even  if  there  is  no  wind.  Consequently,  if  data  for  the  wind  measured 
by  the  pulsed  system  were  plotted  versus  the  actual  wind,  the  plot 
would  not  pass  through  the  origxn.  This  would  represent  a  constant 
error  and  the  system  would  indicate  a  wind  when  in  fact  there  was  no 
wind.  This  has  been  observed  experimentally  at  the  Oregon  Graduate 
Center  test  site. 

The  beam  misalignment  in  the  present  transmitter  can  be  separated 
into  two  parts:  (1)  a  mechanical  slowly  varying  misalignment;  and  (2) 
a  pulse  to  pulse  jitter.  The  mechanical  misalignment  can  be  reduced 
to  an  acceptable  level  on  a  short  term  basis  as  discussed  in  a  section 
on  System  Consideration  by  adjusting  a  steering  mirror  for  one  of  the 
beams  to  make  them  colinear.  Changes  in  ambient  temperature  and/or 
shock  and  vibration  will  bring  the  beams  out  of  alignment.  Typically, 
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the  beams  will  drift  several  spot  diameters  apart  overnight  and  seem  to 
be  in  a  continuous  process  of  drifting  in  and  out  of  alignment  during 
operation  of  the  system. 

The  pulse  to  pulse  jitter  in  the  present  laser  transmitter  is  of 
the  order  of  a  spot  diameter.  It  cannot  be  reduced  with  the  present 
laser  by  adjustment.  Reduction  of  pulse  to  pulse  jitter  will  require 
a  different  approach  to  the  laser  design.  Based  on  the  results  of 
figure  3  and  our  operating  experience,  it  appears  that  a  laser  beam 
alignment  and  RMS  jitter  specification  of  around  one-tenth  the  beam 
divergence  should  allow  good  system  operation.  A  more  definitive 
specification  will  be  available  after  the  analytic  work  on  the  beam 
alignment  problem  is  completed, 

LASER  COHERENCE 

A  spatially  incoherent  but  temporarily  coherent  source  gives  rise 
to  a  vacuum  speckle  of  unity  contrast  ratio.  This  can  represent  a 
strong  source  of  noise  in  such  applications  as  remote  sensing,  optical 
radar,  coherent  imaging  and  COAT  systems  or  it  can  be  the  signal  itself 
as  in  speckle  interferometry.  When  these  applications  involve  propa¬ 
gation  through  the  atmosphere,  the  interaction  of  the  speckle  field 
with  the  turbulent  atmosphere  can  change  its  character.  Consequently, 
an  understanding  of  the  interaction  of  the  speckle  field  with  the  at¬ 
mosphere  is  necessary. 

There  are  nine  combinations  of  spatial  and  temporal  coherence  of 
the  source  as  indicated  in  table  1.  The  strongest  speckle  field,  of 
course,  occurs  for  a  temporarily  coherent  but  spatially  incoherent 
source.  As  the  spatial  coherence  increases  and/or  the  temporal  coher¬ 
ence  decreases,  the  vacuum  contrast  ratio  decreases  from  unity  and 
approaches  zero  as  the  source  either  becomes  spatially  coherent  or 
temporarily  incoherent  as  shown  in  figure  4.  Consequently,  only  four 
of  the  regimes  shov.ui  in  table  1  exhibit  speckle. 

There  are  two  common  physical  configurations  with  regard  to 
speckle  field  propagation  through  the  turbulent  atmosphere.  They  are 
the  double  ended  configuration,  such  as  with  speckle  interferometry, 
where  the  transmitter  is  at  one  end  of  the  propagation  path  and  the 
receiver  is  at  the  othei,  and  the  single  ended  configuration,  which 
applies  to  the  remote  sensing,  optical  radar,  coherent  imaging  and 
COAT  system  application,  where  the  transmitter  and  receiver  are  at  the 
same  end  of  the  path.  In  the  latter  case,  the  effect  of  the  atmos¬ 
phere  on  the  transmitted  energy  must  be  evaluated  both  for  the  propa¬ 
gation  from  the  transmitter  to  the  target  and  for  the  propagation  from 
the  target  back  to  the  receiver.  If  a  temporally  and  spatially  co¬ 
herent  source  is  used,  then  the  speckle  field  source  is  the  energy 
scattered  from  the  target.  Consequently,  the  spatial  coherence  for 
the  return  trip  propagation  is  determined  by  the  target,  with  a  mirror 
or  retroreflector  yielding  a  spatially  coherent  source,  a  perfectly 
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diffuse  target  yielding  a  spatially  incoherent  source  and  with  most 
real  targets  falling  somewhere  in  between  and  giving  both  specular  and 
diffuse  scattering.  The  specular  part  is  usually  referred  to  as  glint 
and  plays  an  important  role  in  COAT  systems. 

Past  work  at  the  Oregon  Graduate  Center  on  speckle  propagation 
through  .he  turbulent  atmosphere  has  been  concentrated  primarily  on 
the  case  of  a  temporally  and  spatially  coherent  transmitter  used  in 
conjunction  with  a  perfectly  diffuse  target.  Analytic  formulations 
have  been  developed  (refs.  17  and  18)  fbr  the  first  and  second  order 
statistics  of  the  intensity.  In  addition,  the  variance  and  covariance 
of  the  intensity  have  been  measured  experimentally  (ref.  19)  and  ex¬ 
cellent  quantitative  agreement  was  obtained  with  the  theory  for  the 
covariance  and  a  good  qualitative  agreement  was  obtained  with  the 
analytic  results  for  the  variance. 

It  was  found  that  three  scale  sizes  are  important  to  the  covar¬ 
iance.  They  are  the  speckle  scale  size  which,  of  course,  depends  on 
the  transmitter  spot  size  on  the  target,  the  atmosphere  phase  per¬ 
turbation  scale  size  p  and  the  atmosphere  log-amplitude  perturbation 
scale  size  /L/k.  At  very  low  turbulence  levels,  the  speckle  scale 
size  is  dominant  and  at  very  high  turbulence  levels  p  is  dominant. 

In  between,  any  of  the  scale  sizes  may  be  dominant  depending  on  their 
relative  size  and  the  integrated  path  turbulence. 

There  is  a  simplified  theory  available  (ref.  17)  for  low  turbu¬ 
lence  that  is  useful  for  some  problems.  It  has  also  been  shown  ana¬ 
lytically  (ref.  17)  and  experimentally  (ref.  19)  that  the  intensity 
statistics  are  exponential  at  low  turbulence  levels. 

The  normalized  variance  for  this  regime  is  unity  at  very  low 
turbulence;  and  as  the  turbulence  increases,  it  rises  above  unity 
until  it  reaches  a  peak  value  of  around  1.25.  As  turbulence  further 
increases,  the  normalized  variance  decreases  until  it  again  reaches 
unity  as  shown  in  figure  3. 

When  the  transmitted  beam  is  scattered  by  the  target,  the  phase 
of  the  scattered  field  is  completely  randomized.  Consequently,  the 
phase  perturbation  caused  by  the  atmosphere  cannot  further  increase 
the  variance.  However,  the  log-amplitude  perturbation  is  not  so  con¬ 
strained  and  accounts  for  the  increase  of  the  normalized  variance 
above  unity.  As  the  turbulence  increases,  the  log-amplitude  pertur¬ 
bation  eventually  saturates  and  the  atmospheric  perturbation  becomes 
dominated  by  the  phase  perturbation.  Consequently,  the  normalized 
variance  returns  to  unity  at  high  turbulence  levels. 

Work  on  the  temporally  coherent,  spatially  incoherent  regime  is 
considered  to  be  substantially  complete.  No  new  work  in  this  regime 
is  being  contemplated  at  this  time  by  the  Oregon  Graduate  Center. 
However,  for  reasons  discussed  below,  of  most  importance  to  the  remote 
crosswind  measurement  system  is  the  temporally  and  spatially,  partially 
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coherent  regime.  Work  Is  in  process  on  this  important  problem. 

For  the  remote  crosswind  measurement  system,  the  fluctuations  in 
the  received  intensity  are  caused  both  by  the  speckle  and  by  the  at¬ 
mospheric  turbulence.  The  speckle  dominates  at  low  turbulence  and  the 
atmosphere  dominates  at  higher  turbulence  levels.  The  transition 
occurs  at  a  turbulence  level  that  decreases  as  the  temporal  coherence 
of  the  laser  decreases. 

The  transition  to  domination  by  speckle  can  be  observed  by  moni¬ 
toring  the  normalized  autocovariance  function.  The  fields  from  the 
transmitter  laser  pulses  are  not  correlated  to  any  great  degree  (if  at 
all)  on  a  pulse  to  pulse  basis.  Consequently,  in  the  absence  of  tur¬ 
bulence,  the  autocovariance  is  near  zero.  As  the  turbulence  increases, 
the  autocovariance  will  increase  due  to  the  common  atmospheric  turbu¬ 
lence  through  which  the  laser  pulses  propagate  until  the  turbulence 
dominates  the  speckle.  This  is  illustrated  qualitatively  in  figure  A 
and  has  been  observed  experimentally  with  the  pulsed  system. 

The  speckle  fields  of  course  carry  no  information  about  the  cross- 
wind.  Consequently,  for  a  given  target  range  there  will  be  a  turbu¬ 
lence  level  below  which  the  system  will  not  respond  to  a  crosswind. 

The  turbulence  level  at  which  this  occurs  decreases  with  decreasing 
laser  coherence.  It  also  is  affected  by  beam  misalignment  (see  section 
on  Laser  Beam  Alignment).  Fortunately,  the  laser  used  in  the  system 
has  a  low  coherence  which  minimizes  the  effect  on  system  operation. 
Also,  at  longer  target  ranges  it  should  be  less  of  a  problem.  It 
should  also  be  noted  that  as  the  speckle  fields  become  dominant,  the 
indicated  crosswind  velocity  will  drop  towards  zero  and  therefore  will 
not  make  the  fire  control  solution  any  worse  than  it  would  be  without 
the  crosswind  correction.  However,  the  transition  point  should  be 
made  as  low  as  possible  to  enhance  the  utility  of  the  system.  Also, 
it  should  be  noted  that  the  autocovariance  can  be  used  as  an  indicator 
of  whether  or  not  the  system  is  accurately  measuring  the  wind. 

The  availability  of  low  turbulence  data  has  allowed  us  to  gain 
some  insight  into  the  relationship  between  the  normalized  variance  and 
the  turbulence  strength  for  the  partially  coherent  transmitter  used  in 
the  system.  The  variance  is  an  important  parameter  in  that  it  is  used 
to  normalize  the  measured  statistics  that  are  used  to  estimate  the 
crosswind.  Normalization  is  required  to  make  the  wind  calibration 
constant  independent  of  turbulence  level. 

The  use  of  the  variance  for  normalization  implies  that  it  is  pro¬ 
portional  to  the  strength  of  turbulence  C  2.  Figure  3  shows  qualita¬ 
tively  the  relationship  between  the  normaSized  variance  and  the 
strength  of  turbulence  for  three  degrees  of  coherence.  The  perfectly 
coherent  source,  diffuse  target  problem  has  been  solved  analytically 
and  the  analysis  has  been  verified  by  experimental  work  as  indicated 
above.  Consequently,  the  relationship  shown  in  figure  3  for  the 


perfectly  coherent  source  is  correct. 

The  partially  coherent  source,  diffuse  target  problem  has  not  been 
solved  analytically;  and  so,  the  relationships  shown  in  figure  3  for 
this  case  are  partly  coniecfure  and  partly  based  on  our  pulsed  experi¬ 
mental  work.  For  the  limiting  case  of  the  incoherent  source,  there  is 
no  residual  variance  due  to  upeckle;  and  so  the  variance  must  be  zero 
in  the  absence  of  turbulence.  In  addition,  there  are  theoretical 
reasons  to  believe  that  for  the  incoherent  case  the  variance  and  the 
turbulence  have  a  linear  relationship  until  saturation  occurs.  Conse¬ 
quently,  the  variance  should  be  proportional  to  the  strength  of  turbu¬ 
lence  for  this  case. 

If  the  laser  source  were  perfectly  incoherent,  then  normalizing  by 
the  variance  would  compensate  for  the  strength  of  turbulence.  However, 
the  laser  source  used  in  the  system  is  partially  coherent  and  conse¬ 
quently,  as  shown  in  figure  3  there  will  be  a  residual  variance  at  zero 
turbulence  due  to  the  speckle  field.  This  does  complicate  the  problem 
somewhat,  but  fortunately  the  relationship  appears  to  still  be  linear. 
However,  our  pulsed  experimental  work  to  date  indicates  that  the  re¬ 
lationship  is  linear  with  respect  to  and  not  C  2. 

In  summary  then,  it  appears  that  the  temporal  coherence  of  the 
laser  transmitter  used  in  the  system  should  be:  (1)  as  low  as  possible 
in  order  to  lower  the  minimum  turbulence  level  at  which  the  system  will 
accurately  measure  the  wind,  and  (?.)  should  remain  constant  in  order  to 
allow  the  system  to  be  normalized  with  respect  to  turbulence  level. 
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Normalized  Autocovariance  nj  Normalized  Variance 


igure  3. 


Relationship  of  the  variance  to  turbulence  strength 
for  different  values  of  laser  temporal  coherence. 


Turbulence  Strength 

Figure  4.  Transition  of  autocovariance  at  low  turbulence  levels. 
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PULSED  EXPERIMENTAL  WORK 


Variance  data  that  we  have  collected  with  a  focused  transmitter 
and  with  a  target  at  500  meters  is  shown  in  figure  5  for  pulse  Hi  and 
in  figure  6  for  pulse  02.  Each  curve  represents  12,200  shots  taken  in 
forty  shot  bursts  (A  seconds).  Data  obtained  for  the  same  turbulence 
level  were  averaged  before  plotting.  However,  the  linear  regression 
analysis  to  obtain  the  best  straight  line  fit  and  the  RMS  error  calcu¬ 
lation  were  performed  on  the  data  before  averaging.  The  data  does  in¬ 
dicate  a  linear  relationship  which  has  very  positive  implications  with 
regard  to  the  design  of  a  pulsed  laser  remote  wind  sensing  system.  It 
also  should  allow  the  system  to  remotely  sense  the  turbulence  level 
which  may  have  some  impact  with  respect  to  characterizing  the  electro¬ 
optic  environment. 

The  linear  fit  to  the  data  yields  almost  identical  equations  for 
each  pulse.  This  indicates  that  at  least  on  the  average,  both  lasers 
in  the  transmitter  had  nearly  the  same  degree  of  coherence.  The  RMS 
errors  for  each  pulse  are  different  however.  This  implies  that  the 
degree  of  coherence  fluctuation  is  different  for  the  two  lasers.  It 
is  important  that  both  lasers  have  nearly  the  same  degree  of  coherence 
and  that  the  degree  of  coherence  is  stable.  This  will  have  to  be  taken 
into  account  in  future  transmitter  designs. 

An  analytic  solution  to  the  partially  coherent  source,  diffuse 
target  problem  is  needed.  It  has  been  proposed  as  a  task  for  the  next 
contract  period. 

Covariance  data  versus  turbulence  strength  for  a  focused  trans¬ 
mitter  with  a  target  at  500  meters  is  shown  in  figures  7  and  8.  A 
theoretical  curve  is  also  shown  based  on  our  earlier  work  on  the  per¬ 
fectly  coherent  source.  Although,  the  actual  source  used  was  only 
partially  coherent,  it  is  expected  that  the  trend  with  increasing 
turbulence  level  should  be  the  same.  As  expected,  the  covariance  both 
for  the  data  and  the  theoretical  curve  falls  off  with  increasing  turbu¬ 
lence  level  indicating  a  decrease  in  scale  size.  The  covariance  does 
have  an  effect  on  both  error  and  calibration  and  is  a  function  of  range 
as  well  as  turbulence  level.  However,  it  can  be  controlled  by  changing 
the  apparent  detector  spacing. 

Crosswind  data  taken  during  the  contract  period  is  shown  in 
figures  9  through  13  as  time  series.  In  addition,  the  data  corres¬ 
ponding  to  the  data  plots  is  given  in  tables  2  through  6.  In  all 
cases,  the  target  is  located  500  meters  from  the  transmitter.  The 
curves  and  data  labelled  Campbell  unit  represent  the  wind  as  measured 
by  a  commercial,  cw  laser  anemometer.  It  is  a  space  averaging  ane¬ 
mometer  and  has  an  accuracy  of  ±  0.5  meter /second.  The  path  weighting 
function  for  the  Campbell  unit  is  not  the  same  as  for  the  pulsed 
system.  Consequently,  part  of  the  difference  in  their  readings  may  be 
attributed  to  this,  which  is  particularly  true  during  low,  changing,  or 
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Standard  Deviation  at  Normalized 
Received  Intensity 


0.4 


Oj  ■  0.142  + .  379  ov, 
RMS  Error  «  13.1  % 


Focused 
300  Mete.'s 
Pulse  #  I 

305  Four  Second  Samples 


0  ■ - ■ - 1 - 1 - 1 _ I 

o  .1  .2  .3  A  .5 

Log -Amplitude  Standard  Deviation 


tigurv.  5.  Standard  deviation  of  normalized  received  intensity 

versus  a  for  Pulse  If  1. 
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Standard  Deviation  at  Normalized 
Received  Intensity 


Normalized  Covariance  of  the  Received  Intensity 


Figure  7.  Covariance,  focused. 
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Normalized  Covariance  of  the  Received  Intensity 


Detector  Spacing  =  4.15  tnra 


Figure  8.  Covariance,  defocused. 
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- Campbell  Unit 

- Pulsed  Unit 

Focused 

Range  =  500  meters 
Cn2  =  2  x  I0"13 
Number  of  Samples  =  40 
t  =  99S  micro  seconds 
4/4/77  Data 
RMS  Difference  =  0.44 


Table  2.  Crosswind  data,  C  2  =  2 

n 

Oregon  Graduate  Center 


Range  =500  meters 

Time  Between  Pulses 

Number  of  Samples  = 

C  2  =  2  x  10-1 3 
n 

Campbell 

0.2 

0.8 

0.9 

0.8 

0.8 

0.9 

1.1 

1.3 

1.4 

1.3 
1.0 
1.0 
0.8 
1.2 
1.2 
0.9 
0.6 
0.1 

0.6 

1  /. 

A*  *T 

1.2 

1.4 

0.6 


=  999  microseconds 
40 


Slope 

Pulsed 

24 

0.39 

-5 

-0.03 

19 

0.31 

2 

0.03 

37 

0.6 

48 

0.77 

68 

1.09 

69 

1.11 

71 

1.14 

67 

1.08 

75 

1.21 

65 

1.05 

54 

0.87 

63 

1.0 

61 

0.98 

101 

1.62 

48 

0.77 

35 

0.56 

28 

0.45 

OO 

4.4. 

r\  i  c 
\j  • 

46 

0.74 

48 

0.77 

15 

0.24 

RMS  Difference  =  .435 
* 

Using  calibration  constant  that  minimizes 
Focused  with  1"  transmitter  stop 


x  10-13. 


Difference 

0.19 

-0.88 

-0.59 

-0.77 

-0.2 

-0.13 

-0.01 

-0.19 

-0.26 

-0.22 

+0.21 

+0.05 

+0.07 

-0.20 

-0.22 

+0.72 

+0.17 

+0.46 

-0.15 

1  AC 

~  JL  • 

-0.46 

-0.63 

-0.36 


the  error  = 


I 


tn 

*o 


puoo9s/sjd*9^  u;  pui^ 


(N 

CM 


II 

eg 


T3 

C 


5b 

•H 

Pn 


K 
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Time  in  Minutes 


Table  3 


.  Crosswind  data,  C^=1.22xl0*3 

n 

Oregon  Graduate  Center 


Range  **  500  meters 

Time  Between  Pulses  =  999  microseconds 

Number  of  Samples  =  AO 

C  2  -  1.22  x  10"13 
n 


Focused  Beam 


pbell 

Slope 

Pulsed* 

1.9 

86 

1.88 

2.1 

120 

2.62 

2.6 

63 

1.38 

2.3 

9A 

2.06 

2.6 

61 

1.33 

1.9 

80 

1.75 

1.5 

5A 

1.18 

1.4 

A1 

0.9 

2.3 

72 

1.6 

2.9 

123 

2.70 

2.7 

119 

2.60 

2.1 

92 

2.01 

1.3 

96 

2.10 

1.8 

110 

2. A0 

1.5 

8A 

1.8A 

1.5 

55 

1.20 

1.3 

38 

0.83 

1.  A 

71 

1.55 

1.5 

83 

1.82 

1.0 

53 

1.16 

1.  A 

89 

1.95 

1.6 

90 

1.97 

RMS  Difference  =  0.54  meters  per  second 
*Using  calibration  constant  =  0.02186 
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Time  in  Minutes 


Table  4.  Crosswind  data,  C  2  =  7,4  x  lO-1*1. 

n 

Oregon  Graduate  Center 

Range  =  500  meters 

Time  Between  Pulses  =  999  microseconds 

Number  of  Samples  =  40 

C  2  =  7.4  x  HT14 
n 


Campbell 

Slope 

Pulsed 

1.8 

114 

1.51 

1.6 

171 

2.26 

1.8 

91 

1.2 

1.6 

95 

1.26 

1.7 

75 

0.99 

1.5 

144 

1.91 

1.5 

105 

1.39 

1.5 

93 

1.23 

1.4 

150 

1.99 

1.9 

57 

0.75 

1.6 

105 

1.39 

1.3 

40 

0.53 

1.5 

130 

1.72 

1.3 

122 

1.62 

W2  =  2.50 
S2  12,553 
WS  =  166 


RMS  Difference  =  0.55 
* 

Using  Calibration  Constant  =  0.01324 
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wind  in  Meters 


3.0  h 


2.0  h 


- Campbell  Unit 

- Pulsed  Unit 

Collimated 

Range  =  500meters 

Cn2  =  6.4  xIO'14 

Number  of  Samples  =40 

r  =  999  micro  seconds 

9/28/77  Data 

RMS  Difference  =0.68 


10  15  20 

Time  in  Minutes 


Figure  12.  Crosswind,  C  2  =  6.4  x  10 
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Table  5.  Crosswind  Data  C  2  =  6.4  x  10-14. 

n 

Oregon  Graduate  Center 

Range  =500  meters 

Time  Between  Pulses  =  999  microseconds 
Number  of  Samples  =  40 
Cn2  =  6.4  x  10"14 
Collimated  Beam 

* 

Campbell  Unit  Slope  Pulsed 


1.7 

62 

0.85 

1.8 

53 

0.72 

1.8 

58 

0.79 

1.2 

70 

0.96 

1.3 

12 

0.16 

0.9 

21 

0.29 

1.2 

51 

0.70 

0.3 

35 

0.48 

0.4 

79 

1.08 

0.3 

60 

0.82 

0 

85 

1.16 

0 

60 

0.82 

0.1 

92 

1.26 

0.4 

38 

0.52 

0.8 

122 

1.67 

1.0 

56 

0.77 

1.4 

108 

1.48 

1.2 

87 

1.19 

2.0 

71 

0.97 

1.9 

155 

2.12 

1.4 

98 

1.34 

1.5 

106 

1.45 

1.8 

42 

0.57 

1.2 

110 

1.5 

1.8 

127 

1.74 

1.8 

132 

1.8 

1.7 

85 

1.16 

1.4 

110 

1.5 

1.3 

132 

1.8 

1.8 

106 

1.45 

1.6 

111 

1.52 

1.9 

79 

1.08 

2.1 

156 

2.13 

2.2 

52 

0.71 

1.7 

133 

1.82 

W2  =  2.039  S2  =  8.45  x  103  WS  =  1.155  x  102 

RMS  Difference  =  0.68  m/second  *Using  Calibration  Constant  =  0.01367 
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Tim*  in  Minuf** 


Table  6.  Crosswind  data,  C  2  ■  3.89  *  10" 1,1 

n 

Oregon  Graduate  Center 


Range  =  500  meters 

Time  Between  Pulses  «*  999  microseconds 
Number  of  Samples  =  40 
C  2  °  3.89  x  10" ^ 

Focused  Beam 


Campbell  unit 

Propellor  anemometer 

Slope 

Pulsed  1 

Pulsed  2 

1.2 

1.9 

63 

1.03 

1.04 

1.6 

1.5 

87 

1.43 

1.43 

1.9 

3.2 

59 

0.97 

0.97 

2.4 

0 

45 

0.74 

0.74 

1.9 

0 

20 

0.33 

0.33 

1.3 

0 

13 

0.21 

0.21 

1.3 

0.3 

35 

0.57 

0.58 

1.1 

0 

36 

0.59 

0.59 

1.6 

0 

69 

1.13 

1.13 

1.5 

0.2 

44 

0.72 

0.72 

1.8 

1.5 

1 

0 

0 

1.6 

1.6 

23 

0.38 

0.38 

1.9 

0.2 

50 

0.82 

0.82 

1.9 

0.2 

16 

0.26 

0.26 

2.1 

1.0 

52 

0.85 

0.85 

1.4 

0.7 

75 

1.23 

1.23 

1.1 

0.1 

57 

0.93 

0.94 

3.3 

1.4 

133 

2.18 

2.19 

0.8 

1.6 

106 

1.74 

1.74 

0.3 

1.8 

132 

2.16 

2.17 

0.9 

1.1 

48 

0.79 

0.79 

0.8 

0 

70 

1.15 

1.15 

0.8 

0 

152 

2. 49 

r%  r  r\ 

i.  .  JU 

1.7 

0.1 

90 

1.4/ 

1.48 

2.6 

5.2 

104 

1.70 

1.71 

2.1 

3.6 

95 

1.56 

1.56 

1.8 

2.8 

111 

1.82 

1.82 

1.5 

1.7 

87 

1.42 

143 

1.3 

0.7 

63 

1.03 

1.04 

1.4 

1.1 

132 

2.16 

2.17 

1.3 

0 

123 

2.02 

2.02 

1.6 

1.4 

178 

2.29 

2.93 

2.3 

1.2 

106 

1.74 

1.74 

1.8 

3.7 

115 

1.88 

1.89 

2.1 

3.5 

103 

1.69 

1.69 

1.8 

3.2 

96 

1.57 

i.58 
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Table 

6.  Crosswind  data,  C 

n 

2  =  3. 

89  x  10-1*4. 

Campbell  unit 

Propellor  anemometer 

Slope 

* 

Pulsed  1 

Pulsed 

2.0 

3.3 

114 

1.87 

1.87 

1.7 

3.2 

96 

1.57 

1.58 

1.6 

2.  7 

106 

1.74 

1.74 

1.7 

0.4 

38 

0.62 

0.62 

1.8 

0.2 

78 

1.28 

1.28 

1.4 

0.4 

14 

0.23 

0.23 

1.3 

0.9 

75 

1.23 

1.23 

* 

Comparison  to  Campbell  Unit 

W2  =  2.857 

S2  =  7.569  x  103 

WS  =  1.24  x  io2 

Calibration  Constant  = 
0.01638 

RMS  Difference  =  0.91 


** 

Comparison  to  Propellor  Anemometer 
W2  »  3.563 
S2  =  7.569  x  ios 
WS  =  1.244  x  io2 
Calibration  Constant  =  0.01644 
RMS  Difference  =  1.232 
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gusty  wind  conditions.  The  column  labeled  ’slope'  in  the  tobies  is  the 
slope  of  the  time  delayed  covariance  function  at  zero  time  delay  and  is 
proportional  to  the  path  averaged  crosswind. 

For  each  data  set,  a  calibration  constant  was  calculated  from  the 
data  such  that  the  mean  square  error  was  a  minimum.  This  constant  was 
then  used  to  calculate  the  pulsed  system  crosswind  and  the  RMS  error 
from  the  slope  and  in  situ  data. 

In  studying  the  results,  it  is  apparent  in  some  of  the  figures, 
particularly  the  first  part  of  figures  12  and  13,  that  there  are 
periods  when  the  pulsed  system  and  the  Campbell  unit  do  not  track  ver> 
well.  During  the3e  periods  th°  system  generally  reads  low.  We  believe 
that  this  is  caused  by  the  two  transmitter  beams  wandering  in  and  out 
of  alignment.  It  shows  up  to  a  lesser  degree  at  higner  turbulence 
levels  which  together  with  the  low  readings  tends  to  support  this 
conclusion. 

During  the  focused  beam  measurements  shown  in  figure  11,  a  single, 
propellor  anemometer  was  used  in  addition  to  the  Campbell  unit.  It 
gives  a  single  point  wind  reading  and  yields  some  insight  into  the  de¬ 
gree  of  fluctuation  of  wind  at  e  single  point  along  the  path.  It  is 
interesting  to  note  that  using  the  single  propellor  a.iemometcr  for  the 
in  situ  data  yielded  the  same  calibration  constant  as  using  the  Camp¬ 
bell  unit  for  the  in  situ  data. 

One  of  the  data  sets  was  taken  with  the  transmitter  beams  colli¬ 
mated  rather  than  focused  (figure  12).  No  substantial  difference  in 
performance  was  noted  between  the  focused  and  collimated  cases.  How¬ 
ever,  there  are  two  reasons  to  believe  '.hat  better  performance,  will  be 
achieved  with  a  focused  transmitter:  (1)  the  low  turbulence  limit  of 
operation  is  determined  by  the  relative  scale  sizes  produced  by  the 
speckle  and  by  the  atmosphere  and  their  relative  scintillation 
strengths.  Since  for  a  focused  transmitter  the  speckle  scale  size  is 
equal  to  “he  transmitter  size,  /b/k,  will  dominate  at  the  ranges  of 
interest  until  the  atmospheric  scintillation  strength  is  reduced  below 
the  speckle  value.  The  speckle  scale  size  from  a  collimated  3-inch 
beam  is  much  smaller  than  the  beam  size;  and,  consequently,  it  is  felt 
that  the  speckle  would  begin  to  dominate  at  a  higher  turbulence  level 
than  with  a  focused  beam.  In  addition,  the  larger  spot  size  on  the 
target  reduces  the  atmospheric  scintillation  strength.  There  may 
however  be  a  trade  off  involving  beam  size  that  wil.  allow  collimated 
beam  performance  to  very  nearly  match  the  focused  beam  performance. 

The  use  of  a  collimated  beam  is  desirable  from  the  standpoint  that  it 
would  simplify  system  design.  (2)  There  is  a  possibility  that  struc¬ 
ture  on  the  target  may  affect  system  performance.  By  using  a  focused 
beam,  the  spot  size  is  determined  primarily  by  the  atmosphere  most  of 
the  time.  Consequently,  any  structure  on  the  target  (within  the  spot) 
will  not  be  observed  by  the  receiver  and  therefore  should  not  have  any 
effect  on  system  performance.  There  may  of  course  be  a  collimated 
beam  size  that  accomplishes  this  same  result. 
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CW  EXPERIMENTAL  WORK 


Cur  cw  experimental  system  (refs.  18,19)  was  used  to  collect  data 
on  the  time  delayed  statistics.  Specifically,  the  autocorrelation 
function  and  the  time  delayed  covariance  runction  were  measured,  and 
good  qualitative  agreement  with  the  theory  was  obtained. 

Measured  autocorrelation  functions  for  two  wind  speeds  are  shown  in 
figures  14  and  15,  The  shape  of  the  curves  is  consistent  with  theory; 
and  as  the  wind  speed  increases,  the  correlation  for  a  given  time  delay 
decreases  as  it  should. 

Figures  16  through  22  show  measured  time  delayed  covariance  func¬ 
tions  for  both  north  and  south  winds.  The  shape  is  consistent  with 
theory,  and  the  peaks  are  on  opposite  sides  of  zero  time  delay  for 
opposing  wind  directions,  which  yields  a  positive  slope  for  a  south 
wind  and  a  negative  slope  for  a  north  wind.  The  slope  data  and  cross- 
wind  calibration  constant  for  each  of  the  curves  is  shown  in  table  7. 

The  measured  slope  at  zero  time  delay  is  reasonably  consistent  with  that 
measured  by  the  pulsed  system. 

Good  quantitative  results  are  difficult  to  obtain  for  the  time  de¬ 
layed  statistics  because  the  wind  is  not  uniform  along  the  path  and 
generally  fluctuates  considerably  with  time.  Consequently,  it  takes 
considerable  time  and  effort  to  obtain  good  data. 

Tn  addition,  more  data  was  collected  for  the  normalized  variance 
of  the  in*-msity.  It  is  shown  in  figure  23  along  with  a  complete 
theoretic; 1  curve  including  the  transition  from  low  to  saturated  turbu- 
.1ence.  The  formulation  for  the  theoretical  curve  was  developing  during 
the  contract  period  and  is  contained  in  the  next  section. 
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AUTOCORR 


TAU 

0.5  millisecond  per  division 


Figure  14.  Autocorrelation  function, 
1.3  raeters/second  crosswind. 
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TAU 

0.5  millisecond  per  division 

Figure  15.  Autocorrelation  function, 

0.7  meter/second  crosswind. 
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Collimated 
Beam  Size  =  1.35  cm 
Wavelength  =  0.488  pm 
500  meters 

Crosswind  =  1.2  m/sec 
Run  #SP4R1 

Detector  Spacing  =  4.5  ram 
C  2  =  3.2  x  10-13 
n 


Figure  16.  Time 
1.2  meters 
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Collimated 
Beam  Size  =  1.35  cm 
Wavelength  =  0. 488  um 
500  meters 

Crosswind  =1.8  m/sec 
Run  //SP2R1 

Detector  Spacing  =  4.5  mm 
C  2  =  6.6  x  10~*4 


TAU 

0.5  millisecond  per  division 

Figure  17.  Time  delayed  covariance  function, 

1.8  meters  per  second  north  wind. 
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Figure  18,  Time  delayed  covariance  function, 
1.2  meters  per  second  south  wind. 


Collimated 
Beam  Size  =  1.35  mm 
Wavelength  =  0.488  um 
500  meters 

Crosswind  =0.7  m/sec 
Run  //SP2R8 

Detector  Spacing  =  4.5  mm 

C  2  -  1.2  x  10~*3 
n 


+  '-'I— . ♦— -»  -  « . ♦  '» 


Figure  20.  Time 
0.7  meter  ; 


0.5  millisecond  per  division 


elayed  covariance  function, 
r  second  north  wind. 
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II. 


Collimated  Beam 
Beam  Size  =  1.35  mm 
Wavelength  =  0.488  pm 
500  meters 

Crosswind  =1.3  m/sec 
Run  // SP2R7 

Detector  Spacing  =  4.5  mm 

C  2  =  1.4  x  10-1 3 
” 


T.\U 

j.5  millisecond  per  division 


Figure  21.  Time  delayed  covariance  function, 
1.3  meters  per  second  north  wind. 
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Table  7.  Slope  data,  t  a  1  millisecond. 


F  or  C 

Run  J? 

Slope 

Crosswind 

Slope/Crossvind  ■=  Cal.  Const 

Collimated 

SP4R1 

+  80 

1.2  S 

0.015 

Collimated 

SP2R1 

-  91.5 

1.8  N 

0.0197 

Focussed 

SP1R2 

62 

1.2  S 

0.019 

Collimated 

SP4R2 

87 

1.4  S 

0.016 

Collimated 

SP2R8 

-  61 

0.7  N 

0.012 

Collimated 

SP2R7 

-  66 

1.3  N 

0.0198 
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Figure  23.  Normalized  variance  of  intensity  versus 


ANALYSIS  AND  NUMERICAL  WORK 


The  theoretical  formulations  (ref  18)  for  the  statistics  of  the 
received  intensity  were  developed  in  terms  of  the  wave  structure  func¬ 
tion  D  and  log-amplitude  covariance  function  C  .  A  direct  calculation 
of  C  ^  for  use  in  the  theoretical  formulations^proved  to  be  too  slow 
to  *  be  useful.  Consequently,  an  algorithm  for  efficiently  generating 
C  has  been  developed.  A  formulation  for  the  log-amplitude  covariance 
function  good  for  all  levels  of  turbulence  is  given  by  (ref.  20) 


cx(p,aT2)  -  2.95  oT2  j  du  [u(l  -  u)J5/& J 
X  exp|-c?T2  [u(l-u)]5/6  f(y)j>  Jq 


dy 


sin2y 

vll/6 


P 

/XL 


— 


where 


a  2  =  .124  k7/6  Ln^6  C  2 
T  n 

03 

f(y)  =  7.02  y5/6  f  dx  x-8/3  [1  -  J  (x)] 
J.7y  ° 


(1) 

(2) 

(3) 


and 


JQ  is  the  zero-order  Bessel  function.  We  first  consider  the  exponen¬ 
tial  factor  f(y),  and  define  the  function  g(q)  as 


g(q)  = J  dx  x-8^3  [1  -  Jo(x)]  dx 

q 

It  can  he  integrated  by  parts,  to  yield 


g(q) 


=  I  q  5^3  U  -  JQ(q)]  +  J  J  dx  x~5/3  J](x) 

q  „ 

J  dx  x“5/3  J  j (x) 


=  f  q_5/3  [i  -  J0(q)3  +  f 


-  j  f  dx 

•'A 


(x) 


(4) 


(5) 


We  can  make  use  of  the  integral  formula 
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and  the  integral  representation 

Jl(bu)  =  j"  dt/l  -  t2  cos  (but) 

T  Jo 

to  reduce  (5)  to 


S(q)  -f  q-5/3  [1  -  J0h»  +  |  2-5/3  iOM 


_6_ 
5  a 


ri  q 

Jdt/l  -  t^  I  dx  x"2/3  cos  xt 

n  n 


o  o 

Now  using  the  integral  identity, 

q 


L 


dx  xP  1  cos  xt  =  ^  [j F !  ( p ;  p  +  1;  iqt)  +  i Fj  (pi ;  p  +  ii 


-  iqt)] 

where  jFj  is  the  confluent  hypergeometric  function  we  obtain 


g(q)  -  T  1 q-5/3  [1  -  j  (q)]  +  1.863893834 


-  i  qi/3  ( dt  /fTli 

n  q  J  L  7  2!  13  4!  19  6! 


+  la 


8h8 


25  8! 


The  integration  over  t  can  now  be  evaluated  by  making  use  of  the 
integral  identity 


/d 

*'  n 


dt  t 


2n 


A~^T2  =  i 


1  r(~2t"1)r(3/;> 


T(n  +  2) 


Performing  the  integration  Eq.  (8)  beet 

,1/3  q7/3  q 1 3/3  q 1 9/3  q25/3  | 

/  /  /  n  **”  ennr  n  e  oni  /  r  /  n  r  ft  /  /  AAAA  1 


g 


(q)  -  3  | 


.37278 


448  29952  2802664  368640000) 

+  A(q) 


where  |  A (q )  |  <  .00041  for  q  <_  —  a. 

For  large  values  of  q  we  can  obtain  an  asymptotic  series  by 

the  asymptotic  expansion  for  J  sx)  in  Eq.(4).  It  becomes 

o 


(6) 

(7) 

(8) 

(9) 

(10) 

using 
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00 


g(q)  ~  y  q“5/3 


-  —  /  dx  x~19/f  cos  x 

q 


03 

-  —  /  dx  x 


“19/6  s£n  x 


(ID 


Using  the  integral  identity 

J  sin  (ax)  =  [iF^p;  \i  +  1;  ia)  -  iFj  (y ;  p  +  1;  -  ia)]  (12) 


Eq. (11)  becomes 


g(q)  ~  f  q  5^3  — — 
2/ii 


'  13tt 
e  12  r 


13tt 

12 


1 

2^ 


'  .  13ir 

e  ^  r 


.13tt 

("  T’  iq)  "  e  12  r(-T-  -  ii). 


(13) 


where  r(a,b)  is  the  incomplete  gamma  function.  Equation  (13)  can  be 
rewritten  as 


,  ^  _  3  -*/,  (1  +  i)  12 

g(q)  ~  t  q  - ~r~  e 

5  2^ 


.  13tt  -  iq 

^  i  rt 


„{19  19  .  \ 

ulT’  T’ iq/ 


.  13tt  . 

1  -  i  12  +  iq 
-  e 

2/n" 


19  19  .  \ 

uIt>  t-  - iq/ 


|  q‘5/3  -jZ  cos  (<]  +  i)  [q  19/6  -  ^  '  T  ^ 

'] 

V2  .  j  7T \  Pl9  —26/6  19  25  31  —  37/g 

-sin(q+I)  [Tq  25/6  _q  37/6 


+  Ii  .  25  .  31  .  37  -43/6  _ 
6  6  6  6  q 


+  ii  .  25  .  31  .  37  .  43  -49/6  _ 


•] 


(14) 
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where  U  is  confluent  hypergeometric  function  of  the  second  kind. 

Using  Eq.(14)  in  (13)  yields 

f(y)  -  9.45  (.7y)5/6  g(.7y)  (15) 

By  using  ( 1.5)  in  the  expression  for  C  given  by  Eq.(l)  and  performing 
the  u  and  y  integrations  using  Gaussian  quadrature,  numerical  values 

of  C  can  be  efficiently  obtained. 

X 

In  addition  to  t^g  above  work,  some  additional  analytic  work  was 
accomplished  on  the  n  order  mutual  coherence  function  first  formu¬ 
lated  in  reference  18.  The  work  has  been  published  (ref.  21);  and  so, 
it  will  not  be  repeated  in  this  report. 
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